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■ Abstract 

^ ■ We analyze the Scarf potential, which exhibits both discrete energy 

bound states and energy bands, through the quantum Hamilton-Jacobi 
approach. The singularity structure and the boundary conditions in the 
CM ■ above approach, naturally isolate the bound and periodic states, once 

the problem is mapped to the zero energy sector of another quasi-exactly 
solvable quantum problem. The energy eigenvalues are obtained without 
having to solve for the corresponding eigenfunctions explicitly. We also 
demonstrate how to find the eigenfunctions through this method. 

CM ; 

q ; 1 Introduction 

Study of periodic potentials has evoked renewed interest in the literature in light 
"^h , of their appearance in Bose-Einstein condensates pQ, (2> and photonic crystals 

Q H 1 [3], @|, 0, [Sj. It has been possible to experimentally change structure of the 

| ' , , potential, so as to produce superfluid - insulator transition 0, the former hav- 

£^ ' ing delocalized states and the latter localized ones. In this context, traditional 

Kronig-Penny model 8 is used for illustrative purposes, wherein the known 
wave functions lead to transcendental equations involving energy and momen- 
tum, when appropriate boundary conditions are implemented. The possibility 
of investigating superfluid-insulator type transitions mentioned above does not 
arise here, due to the lack of any control parameter. Quite sometime back, Scarf 
jj] ■ showed that a solvable model exists, which exhibits both discrete bound states 

and band spectra UJ, as a function of the coupling parameter. The group the- 
oretical aspects of this problem have recently been investigated [TJ]|. The fact 
that Scarf potential yields both bound states and band structure, as a function 
of a coupling parameter, makes this model an ideal one to study the interplay 
of these two types of distinct behavior in a given quantal problem. 

The goal of this paper is to first map the Scarf eigenvalue problem into 
the zero energy sector of another quasi-exactly solvable (QES) problem. We 
then use the quantum Hamilton-Jacobi (QHJ) approach which naturally 
takes advantage of the singularities of the new potential, to isolate the domains 
corresponding to discrete and band spectra. The subtle aspects of the boundary 
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conditions in quantum mechanics, which lead to the existence of both bound 
states and band structure in the Scarf potential, come out naturally in this 
approach. We then proceed to obtain the eigenvalues and eigenfunctions, for 
both the cases. In this procedure, the energy eigenvalues can be obtained, 
without finding the eigenfunctions. The QHJ formalism, being formulated in the 
complex domain where the non-linear Riccati equation replaces the Schrodinger 
equation, makes use of powerful theorems in complex variable theory to obtain 
the solutions. 

Apart from the fact that QHJ formalism is relatively new and requires de- 
tailed study, this approach may provide a different perturbative treatment for 
the traditional problems. As will be clear form the text, WKB approximation 
scheme is close to this method ^2] . In the following section, we briefly describe 
the working principles of the QHJ formalism, which is then used for the anal- 
ysis of the Scarf potential in Section 3. The origin of the bound and the band 
spectra is then illustrated, without getting into the explicit computation of the 
eigenvalues, whose details are given in Section 3.2. We obtain the solutions 
pertaining to both the spectra. We conclude in the final section after pointing 
out various directions for future investigations. 



2 Quantum Hamilton - Jacobi formalism 

The QHJ formalism, formulated as a theory analogous to the classical canonical 
transformation theory ^3], JSJi was proposed by Leacock and Padgett 
in 1983. It has been applied to one dimensional bound state problems and 
separable problems in higher dimensions In our earlier studies, we have 

shown that one could use the QHJ formalism to analyze one dimensional exactly 
solvable (ES), quasi - exactly solvable models, consisting of both periodic and 
aperiodic potentials [H], [T7j, [H], EH, EDI, E3 and the recently discovered PT 
symmetric potentials [22] • The advantage of this method lies in the fact that 
it requires a modest understanding of basic quantum mechanics and complex 
analysis as a prerequisite. 

In this formalism, the logarithmic derivative of the wave function ip{x), given 

by 

p = — ih— \n.tp{x), (1) 
ax 

plays an important role. This is referred to as the quantum momentum function 
(QMF), since it is defined analogous to the classical momentum function as, 
p = Here, S is the Hamilton's characteristic function which is related to 
the wave function by ip{ x ) = &xp(iS/K). Substituting ip(x) in terms of S in the 
Schrodinger equation, 

_ ^^ + VW(l)= ^ (l) ' (2) 
and using the relation between p and ifi, one obtains the non-linear Riccati 
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equation: 




2m{E-V{x)). 



(3) 



The above equation is known as the QHJ equation; here x is treated as a complex 
variable, thereby extending the definition of p to the complex plane. We show 
that one can arrive at the required results by studying the singularity structure 
of the QMF. 
Singularity structure 

The QMF has two types of singularities, the moving and the fixed singularities. 
From one can see that, the n nodes of the n th excited state, whose locations 
depend on the initial conditions and energy, correspond to the singularities of 
p. These are known as the moving singularities. It is a fact that only poles can 
appear as moving singularities in the solutions of the Riccati equation. One can 
calculate the residue at a moving pole xo, where V(x) is analytic, by doing a 
Laurent expansion of p around Xo as, 



Substituting this in J3J and comparing individually the coefficients of different 
powers of x — xq, one obtains 1 = 1, with the corresponding residue equalling 



The fixed singularities originate from the potential and are present in all the 
solutions of the Riccati equation. One can calculate the residue at the fixed 
poles in the same way, as is done for the moving poles. Owing to the quadratic 
nature of the QHJ equation, one obtains two solutions. In order to arrive at the 
right solution, one needs to choose the residue that gives the correct physical 
behavior. The right value of the residue is chosen by applying the appropriate 
boundary conditions, details of which will be given in the text, as and when 
required. Thus, knowing the singularity structure of QMF and the behavior of 
p at infinity, one gets the complete form of the QMF. In all the models studied 
so far, including the periodic potentials, the assumption that, the QMF has 
finite number of singularities, is equivalent to saying that the point at infinity is 
an isolated singular point, has been found to be true. We expect it to be valid 
for the present case also. 

For most exactly solvable models, the QMF has been found to be a rational 
function. As is known, for a rational function the sum of all residues including 
that at infinity is zero. This result has been used to obtain the energy eigenvalues 
for all the models studied in the QHJ approach. It should be pointed out that, 
this condition is equivalent to the exact quantization condition satisfied by the 
action J [TT] : 



Jc 

Hence, for the case of Scarf potential, one first tries to bring the QMF into a 
rational form through a suitable change of variable, as discussed in the next 
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Figure 1: Scarf potential, with s = 2, a = 1, allows bound states. 



section. It is interesting to note that in the classical limit, 

p^p c = y/2m(E-V(x)), 



(6) 



where p c is the classical momentum. The QHJ quantization condition then leads 
to the WKB approximation scheme. The boundary condition @ was originally 
used by Leacock and Padgett to obtain the constraints on the residues [TT]. 



3 The Scarf Potential 

The Scarf potential is given by 

where, a is the potential period. One finds that in the range s > 1/2, the 
potential is an array of infinite potential wells as shown in Fig 1. 

A quantum particle is then confined to only one well, implying that the 
wave function should vanish at x = ±a. Thus, in the above range, the potential 
exhibits bound state spectra. As shown in Fig. 2, in the range < s < 1/2, the 
potential is similar to that of a potential in a crystal lattice, 

leading to the possibility of energy bands. In this scenario, a particle can 
escape to infinity. Therefore the wave function need not vanish at x = ±a. 
However, ip(x) should not diverge anywhere, on physical grounds. 

The QHJ equation for the Scarf potential, with p — —iq and % = 1 in (JSJ, is 
given by, 



n 2 




where A 2 = 2mEa 2 /it 2 . We perform a change of variable using 

y = cot (— ) , (9) 



4 




Figure 2: Scarf potential, with s = 0.4, a = 1, allows band structure. 



which transforms JEt to 



^/i , 2, (/ '/ ~~ /',2 , 1 2,,-, 2 



- -a + r)^ + ^ ( A + (i - s )(! +v))= °- ( 10 ) 

In order to get all the coefficients in the above equation to a rational form, which 
in turn will easily yield the singularities of the QMF, we use the transformation 
equations 

«? = -— (l + y 2 )0; ^x--* (11) 
a 1 + y 

This leads to the QHJ equation in terms of \ as, 

X 2 + ^ + ^— ^ + ^^ = (12) 

X + dj, + (J/2 + 1)2 + + ! ^ 

Henceforth, the above equation will be treated as the QHJ equation and x 
as the QMF. It is interesting to note that, substituting \ — j'Q&i'&iy)) m 
the above equation, one gets a Schrodinger equation, which describes the zero 
energy sector of the potential, (A 2 — l)/(y 2 + l) 2 + (j — s 2 )/(y 2 + l). By analyzing 
the singularity structure of this quasi-exactly solvable problem, we obtain the 
required results for the solvable Scarf potential, as shown below. 

3.1 Form of the QMF x 

The QMF has n moving poles with residue one on the real line, as is clear from 
the Riccati equation. From 1)12(1. one can see that \ has fixed poles at y = ±i. 
Making use of the assumption that the QMF has finite number of moving poles, 
one can write x i n the rational form, separating its analytical and singular parts 
as 

x = ^- + ^.+± — + Q- (13) 
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Here, b\ and b\ are the residues at y = i and y — —i respectively and the 
summation term describes the sum of all the singular parts coming from the 

moving poles. Note that, Y^k=o y -y k = p" faj » wnere Pn(y) is an n th degree 
polynomial. The quantity Q represents the analytic part of x and from (|12[1 one 
can see that \ is bounded for large y. Thus, from Liouville's theorem, Q is a 
constant; denoting it as C, l|13ll can be written as, 

One can calculate the residues at the fixed poles y = ±i, by making a Laurent 
expansion of x around the pole. For example, to calculate the residue at y = i, 
we expand x as, 

X = — — + a + oi (y - i) H . (15) 

2/ - * 

Comparing the coefficients of different powers of (y — i) individually, one 
obtains 

Similarly the other residue at y — —i is found to be 

K - i±±. (17) 

To find the eigenvalues, we now make use of the fact that, for a rational function, 
the sum of all the residues equals zero. As noted earlier, this is equivalent to 
the quantization condition (0 of Leacock and Padgett. Thus, we obtain 

b 1 + b[ + n = d 1 , (18) 

where d\ is the residue at infinity, which is calculated by taking Laurent expan- 
sion of x around the point at infinity: 

X = d + — + ^ + --- . 19 

y y 

Substitution of the above in the QHJ equation yields, 

df - di + (| - s 2 ) - 0, (20) 

from which, the values of di can be deduced: 

l±2s 

dl = — J— ■ (21) 
Substituting the values of the residues in (fTH|l . one obtains 

n=~±sT\ (22) 
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which gives the degree of the polynomial P n {y) in 1141) . From the definition of 
A, one can see that if E < 0, A becomes imaginary, in which case (|22() will not be 
satisfied. Thus, from the above equation, we have the condition E > 0, which 
in turn implies A > and real. Hence, for any range of s, the energy eigenvalues 
are greater than zero. With this condition on A, we now proceed to select the 
values of the residues at the fixed poles and at infinity, which will give us the 
physically acceptable results. 



3.2 Choice of the residues 

One needs to use the boundary conditions obeyed by the QMF to choose 
the right value of residues. Although, there are several ways of implementing 
the boundary conditions in the QHJ formalism, we have chosen the one closest 
to the conventional approach for clarity. First, we shall fix the value of the 
residue at infinity. From the prior discussion of the potential, we know that the 
wave functions should not become infinite anywhere, in particular, for x = ±a. 
From QJ, one can obtain 4>( x ) m terms of the QMF. Writings = — iq and doing 
the change of variable, one obtains the wave function: 



a 

7r y 1 + y 1 



m = exp ( - / - ( ^3 )dy). (23) 



Using the transformation equations in (|llf) . the above expression for the wave 
function becomes 

V>(y)=ex P (y (x-Y^dvy (24) 



M\ua\jui i >i \ < )i u ami '< i a^ \ "~ ' ! 

substituted in (|24|l . yields, 



For large y, the leading behavior of x is obtained as \ ~ which when 



'•<//.)— -M>( / (7-1^7) ''■'') l2r » 



y dl 



( y a + 1 )i/a- 

Using the value of d\ from 121() in the above equation, one obtains 



(26) 



For < s < i, one can see that 0, in the limit y — * oo, x — > ma, with m 
being an integer, for both the values of d\. This range corresponds to the case 
where the potential exhibits band structure. 

For s > 1/2, ip(y) 0j m the limit y — > oo,x — > ma, with m being an 
integer, only if d\ takes the value 1/2 — s. In this way, the two different ranges 
of the potential parameter s emerge simultaneously, while fixing the values of 
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d\ . In order to select the values of b\ and b[ , we note that the bound state and 
band edge wave functions of one dimensional potentials are non-degenerate and 
have definite parity. Parity operation requires that x(~y) — ~x(y)j which in 
turn gives 

h=b[. (28) 
With the above constraint, H18|l becomes 

26i + n = di. (29) 

Finiteness of the wave function as i — > oc, gives the values of d\ in the two 
ranges as 

( l±2s for < s < 1/2, 
di = { ' (30) 

[ i=ge- for a>l/2. 

From the parity constraint, one obtains the restriction on the values of the 
residues at the fixed poles as b\ = b[. Using these results, we proceed proceed 
to calculate the solutions for the two ranges. 



3.3 Case 1 : Band spectrum 

In the range < s < 1/2, we have seen that d\ can take both the values of the 
residues. Taking all the possible combinations of the residues, with b\ = b[ and 
substituting them in (|18J) . we evaluate n, the degree of the polynomial P n (y). 
There are four combinations forming four different sets, as given in the fifth 
column of table I. Since n needs to be positive, from table I, we pick only those 
sets which give a positive integral value for n. As seen earlier, A is a positive 
real constant. Thus, only the sets 1 and 2 will yield positive values for n and 
hence; the other two sets are ruled out. Taking the values of b\ and d\ from 
the sets 1 and 2, substituting them in H29|) and using the definition of A and s, 
we obtain the expressions for the energy eigenvalues corresponding to the two 
band edges of the n th band as, 

Here, correspond to the upper and lower band energies of the n th band. 
These results match with the solutions given in |5] and ^H] . 
The corresponding wave functions follow from 114(1 and 12411 : 

^(y) = {y 2 + \) b ^P n {y). (32) 

To obtain the expression for the polynomial, we substitute \ from ITljl in the 
QHJ equation, which gives a second order differential equation: 

?n(y) + (^) K(v) + 

/1/4-s 2 4fe^ 2 + 26i(l-j/ 2 ) + A 2 -l 

I y 2 + i + (y 2 + 1) 2 



J P n (y) - 0. (33) 
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The sets 1 and 2 have b\ = (1 — A)/2, which yields 

™ + ( wO ™ + A pm(y) = °- (34) 

From Q31[l . one can see that A has two values A = n ± s + |. Substituting these 
in the above equation, one obtains two differential equations corresponding to 
the two energy eigenvalues : 

(y 2 + l)P^(y) + (1 - 2n T 2s)vK(v) + «(« ± 2s)P„(y) = 0. (35) 

Defining y = it, the above equation takes the form of the well known Jacobi 
differential equation 

{l-t 2 )P l ;{t) + {v 1 -u 2 -t{v l + v 2 + 2))P l n {t) + n{n + v 1 + v 2 + l)P n (t) = 0, (36) 

with V\ = v 2 = — n^fs — 1/2, for the corresponding two A values. The expression 
for the two band edge wave functions for the n th band are given by, 

^{y) = (y 2 + l)-^P^(-iy), (37) 
with their respective v\, v 2 values corresponding toA = n±s + l/2. 

3.4 Case 2 : Bound state spectrum 

We proceed in the same way as in case 1 i.e., take all possible combinations of 
bi, b[ and d±, keeping bi = b[ in l|29|l . Since d\ can take only one value 1/2 — ,s, 
only two sets are possible here. Out of these, the set corresponding to b\ = b[ = 
(1 — A)/2 alone, will give a positive value for n. Thus, substituting these values 
of residues in 129|) , one obtains the expression for the energy eigenvalue as 

71-2 ( 1 A 2mV a?\ 

^ = w(2 +n + Vl-^J ' (38) 

where n can take positive integral values. Proceeding as above one obtains the 
Jacobi differential equation in terms of t for the polynomial part : 

(1 - t 2 )P^{t) - 2t(-n - s + l)P;(t) - n(n + 2s)P n (t) = 0. (39) 
The expression for the wave function is then given by 

V>(y) = 0/ 2 + ir^ S2 Hy), (40) 

where s± = s 2 = —n — s — 1/2. 

Hence, as pointed out in the beginning, the two different sectors of the 
Hamiltonian, as a function of the coupling parameter and the eigenvalues emerge 
from general principles of QHJ formalism, relying on the singularity structure of 
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the QMF function. The wave functions corresponding to the definite eigenvalues 
are obtained at the end, which match with the known results |23|. 
Conclusions 

We have mapped the entire Scarf problem, containing the bound state and 
energy bands, to the zero energy sector of a different Hamiltonian, which is 
quasi-exactly solvable. This was achieved through point canonical transforma- 
tions which led to the redistribution of singularities in the complex domain. The 
singularity structure of this new Hamiltonian is transparent enough to clearly 
isolate two different regimes, as a function of the coupling constant. When re- 
lated to the original problem, they turn out to represent discrete levels and the 
band edges. It will be interesting to carefully analyze the equilibrium structure 
of the classical electrostatics problem, associated with the QES system, which 
leads to both bound states and band structure in the quantum domain. In light 
of the current interest in periodic potentials in BEC and photonic crystals, we 
hope the quantum Hamilton- Jacobi based treatment presented here is not only 
illuminating, but may also lead to development of new perturbative treatments 
for non-exactly solvable problems. 
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Table 1. All possible combinations of residues in the range 

< s < 1/2. 



Set 


bi 


b[ 




n = di — 6i — 6j 


remark 


1 


l-A 
2 


l-A 
2 


di 


= 1/2-s 


A-s-i 


A > (s + i) 


2 


l-A 
2 


l-A 
2 


di 


= 1/2 + 8 


A + s - I 


A>-(s-i) 


3 


1+A 
2 


1 + A 
2 


di 


= 1/2-s 


-A-s-i 


not valid 


4 


1+A 
2 


1 + A 
2 


di 


= 1/2 + a 


-A+s-i 


not valid 
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